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Abstract. Given a Poisson process on a bounded interval, its random geo- 
metric graph is the graph whose vertices are the points of the Poisson process 
and edges exist between two points if and only if their distance is less than a 
fixed given threshold. We compute explicitly the distribution of the number 
of connected components of this graph. The proof relies on inverting some 
Laplace transforms. 



1. Motivation 

As technology goes on [TJ [5J|3], one can expect a wide expansion of the so-called 
sensor networks. Such networks represent the next evolutionary step in building, 
utilities, industrial, home, agriculture, defense and many other contexts 0]. 

These networks are built upon a multitude of small and cheap sensors which 
are devices with limited transmission capabilities. Each sensor monitors a region 
around itself by measuring some environmental quantities (e.g., temperature, hu- 
midity), detecting intrusion, etc, and broadcasts its collected informations to other 
sensors or to a central node. The question of whether information can be shared 
among the whole network is then of crucial importance. 

Many researches have recently been dedicated to this problem considering a 
variety of situations. It is possible to categorize three main scenarios: those where 
it is possible to choose the position of each sensor, those where sensors are arbitrarily 
deployed in the target region with the control of a central station and those where 
the sensor locations are random in a decentralized system. 

The problem of the first scenario is that, in many cases, placing the sensors is 
impossible or implies a high cost. Sometimes this impossibility comes from the fact 
that the cost of placing each sensor is too large and sometimes the network has an 
inherent random behavior (like in the ad hoc case, where users move). In addition, 
this policy cannot take into account the configuration of the network in the case of 
failure of some sensor. 

The drawback of the second scenario is a higher unity cost of sensors, since each 
one has to communicate with the central station. Besides, the central station itself 
increases the cost of the whole system. Moreover, if sensors are supposed to know 
their positions, an absolute positioning system has to be included in each sensor, 
making their hardware even more complex and then more expensive. 

It is thus important to investigate the third scenario: randomly located sensors, 
no central station. Actually, if we can predict some characteristics of the topology 
of a random network, the number of sensors (or, as well, the power supply of 
them) can be a priori determined such that a given network may operate with high 
probability. For instance, we can choose the mean number of sensors such that, if 
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they are randomly deployed, there is more than 99% of probability the network to 
be completely connected. 

Usually, sensors are deployed in the plane or in the ambient space, thus math- 
ematically speaking, one has to deal with configurations in R 2 , R 3 or a manifold. 
The recent works of Ghrist and his collaborators [SJ [5] show how, in any dimension, 
algebraic topology can be used to compute the coverage of a given configuration of 
sensors. Trying to pursue their work for random settings, we quickly realized that 
the dimension of the ambient space played a key role. We then first began by the 
analysis of dimension 1, which appeared as the most simple situation. There is here 
no need of the sophisticated tools of algebraic topology. However, it doesn't seem 
that the problem of coverage on a finite length interval has already been solved in 
the full extent we do here. Higher dimensions will be the object of forthcoming 
papers. 

We here address the situation where the radio communications are sufficently 
polarized so that we can consider we have some privileged dimension. Random 
coverage in one dimension has been already studied in different contexts. Some 
years ago already, several analysis were done on the circle ([7J [5] and references 
therein) for a fixed number of points and uniform distribution of points over the 
circle. The question addressed was that of full coverage. More recently, in [5], 
efficient algorithms to determine whether a region is covered considering the sensors 
are deployed over a circle and distributed as a Poisson point process are given. In 
[TUTU], the distribution of a fixed number of clusters (see below for the definition) 
is given. In sensors are actually placed in a plan, have a fixed radius of 

observation. The trace of the covered regions over a line is then studied. 

Our main result is the distribution of the number of connected components 
for a Poisson distribution of sensors in a bounded interval. Our method is very 
much related to queueing theory. Indeed, clusters, i.e., sequence of neighboring 
sensors, are the strict analogous of busy periods. As will appear below, our analysis 
turns down to be that of an M/D/l/1 queue with preemption: when a customer 
arrives during a service, it preempts the server and, since there is no buffer, the 
customer who was in service is removed from the queuing system. To the best of 
our knowledge, such a system has never been studied but the usual methods of 
Laplace transform, renewal processes, work perfectly and with a bit of calculus, 
one can compute all the characteristics we are interested in. 

The paper is organized in the following way: Section II presents the physical 
and random assumptions and defines the relevant quantities to be calculated. The 
calculations and analytical results are presented in Section III. In section IV, two 
other scenarios are presented, considering the number of incomplete clusters and 
clusters placed in a circle. In Section V, numerical examples are presented and 
analyzed. 

2. Problem Formulation 

Let L > be the length of the domain in which sensors are located. We assume 
that sensors are distributed according to a Poisson process of intensity A. Let 
(Xi, 1 < i < n) be the positions of the sensors. We thus know that the random 
variables, AX; = Xj+i — Xi are i.i.d. and exponentially distributed. Due to their 
technological limitations, each sensor can communicate only with other sensors 
within a range e: two sensors, located respectively at x and y, are said to be 
directly connected whenever \x — y\ < e. For i < j. two sensors located at Xj 
and Xj are indirectly connected if Xi and Xj+i are directly connected for any 
I = i, ■ ■ ■ , j — l.A set of sensors directly or indirectly connected is called a cluster 
and the connectivity of the whole network is measured by the number of clusters. 
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The number of points in the interval [0,ir] is denoted by N x = X)^°=o 1{x„<i}- 
The random variable Ai given by 



A, 



Xi if i = 1, 

va£j{Xj\Xj > Ai-i,Xj - Xj-i > e] if i > 1, 



represents the beginning of the i-th cluster, denoted by Ci. In the same way, the 
end of this same cluster, Ei, is defined by 

E t . = infj{X; + e\Xj > A t ,X J+1 - Xj > e}. 

So, the i-th cluster, Ci, has a number of points given by Ne ( — N~Af We define 
the length Bi of Cj as Ei — Ai. The intercluster size, Di, is the distance between 
the end of Cj and the beginning of Cj_|_i, which means that Di = Ai + \ — Ei and 
AAi is the distance between the first points of two consecutive clusters Ci, given 
by AA4 = A i+1 - At = B. L + D. t . 

Remark 1. With this set of assumptions and definitions, we can see our problem 
actually as an M/D/l/1 preemptive queue, Fig. [H In this non- conservative system, 
the service time is deterministic and given by e. When a customer arrives during 
a service, the served customer is removed from the system and replaced by the 
arriving customer. Within this framework, a cluster corresponds to what is called 
a busy period, the intercluster size is the idle time and Ai + Di is the length of the 
i-th cycle. 
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FIGURE 1. Queueing representation of the proposed problem. A 
down arrow denotes that user i starts to be served. An up arrow 
indicates that user i leaves the system without have finished the 
service. A double up arrow illustrates that the service of user % 
finishes. It is also shown the beginning and the end of the ith busy 
period, respectively, Ai and E%. 



The number of complete clusters in [0, L] corresponds to the number of connected 
components f3o(L) (since in dimension 1, it coincides with the Euler characteristics 
of the union of intervals, see [T^l [S]) of the network. The distance between the 
beginning of the first cluster and the beginning of the (i + l)-th one is defined 
as Ui = A^fc- We also define AA"o = Dq = X\. Fig. [2] illustrates these 

definitions. 

Lemma 1. For any i G N!J_ ; Ai and Ei are stopping times. 



4 L. DECREUSEFOND AND E. FERRAZ 




U 2 



Figure 2. Definitions of the relevant quantities of the network: 
distance between points, distance between clusters, the size of clus- 
ters, the size interclusters, the beginning of clusters and the end of 
clusters. 



Proof. Let us consider the filtration T t = o-{N a , a < t}. For i = 1, we have 

{A l <t}^ {X 1 <t}^ {N t > 1} e Tt- 
Thus, Ai is a stopping time. For A 2 , we have 



{A 2 > t} & |J lN t = n, |J I AX, > e, (J {AX k < e} \ \ e T t , 

n>l { j=l { k=j+l 

so A 2 is also a stopping time. We proceed along the same line for others Ai and as 
well for Ei to prove that they are stopping times. □ 



Since N is a strong Markov process, the next corollary is immediate. 

Corollary 2. The set {E>i,Di \ i => 1} is a set of independent random variables. 
Moreover, Di is distributed as an exponential random variable with mean 1/ A and 
the random variables {Bi \ i > 1} are i.i.d. 



3. Calculations 

Theorem 3. The Laplace transform of the distribution of Bi, is given by 

E[e-'*'] = A + * ■ 
Proof. Since AXj is an exponentially distributed random variable, 

E [e- sAX n {AXj < e} ] = £ e-°*\e-*dt = -A_ (l - e -(-+A)«) . 



ON THE ONE DIMENSIONAL POISSON RANDOM GEOMETRIC GRAPH 



Hence, the Laplace transform of the distribution of B\ is given by 

oo 

E [e- sBl ] = E [e~ sBl , N El = n] 



oo 

E ] 

n=l 

■DC 



n-1 



\^ J = i 3 t a {AX „> £ } _[_[ 1{AX,<, 



3=1 



£ (E [e^l^])"- 1 ! [e- sA ^l {A x„> 6} ] e" 
A 



n=l 

oo 

E 

n=0 



s + A 



(l_ e -(«+^)) e" sA e 



-S\ _ — S€ 



X + s 



se \e e se + y 

Using Colloraryj^l we have E [e~ sSl ] = E [e~ sSi ] , which concludes the proof. □ 

From this result, we can immediately calculate the Laplace transform of the 
distribution of AA,. Since AA t = B t + D h we have E [e~ sAA '] = E [ e - s (- B *+- D 0] 
and using Corollary [5] 



E [e- sAAi ] = E [e~ sB *] E [ e -" Di ] 



A 



A + se( x+s ) £ 



Corollary 4. The Laplace transform of the distribution of U n , for n > is given 
by 

X n 



E e 



-sU„ 



(A + se( A + s ) £ )' 



Proof. We use Corollaries [5] and [3] to calculate the Laplace transform of the distri- 
bution of U n , since U n = Y,i=i( B i + A): 

n 

E[e~ sU "] = l[E[e- sB *]E[e- sD <] 



X + s 



X 



a 



A + se( A+s ) c / + 

hence the result. 

Let us define the function p n as 

p n : x 6 R + i — > p n (x) = Pr(/3 (x) = n), 

i.e., p n (x) is the probability of having n clusters in the interval [0, x]. Since for all 
x S R+, < p n (x) < 1, the Laplace transform of p n with respect to ir, 

y* oo 

£{p™}( s ) = / e~ sx p„(2;) da;, 
Jo 

is well defined. 

Theorem 5. for any n > 0, i/ie Laplace transform of p n is given by 

A n e (A+ S ) £ 



(1) 



£{Pn}(s) 



( se (A+ S ) £ + A )' 
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Figure 3. Illustration of the condition equivalent to /3q > n. 



Proof. We note that, see Figure [31 

f {AX + U n -i + B n < L} if n > 1, 
|{AXo<oo} ifn = 0. 



{Po(x) > n} 
Hence 



Pr{p {x) = 0) = 1 - Pr(AX + B 1 <x), 

and 

(2) Pr(A,(a:) = n) = Pr(AX + t/„_i + B n < x) - Pr(AX + f/„ + S n+1 < a:). 
Let 

AX + U n -i+B n ifn>l 
if n = ' 



y„ = 

then we have: 

£{Pr(F n < .)}(«) 



/•OO 

/ Pr(F„ < x)e~ sx dx 
Jo 

/>oo />£ 

/ / dP Yn {y)e~ sx dx 
Jo Jo 



= -E [e- sY «l 

(3) = I A " 

for n > 1, where we used Corollary [5] in the third line. For n — 0, the Laplace 
transform is trivial and given by £{Pr(Y~ < 0}( s ) = V s - Substituting Eq. ([3]) in 
the Laplace transform of both sides of Eq. © yields: 

£{p n }{s) = C{Pv(Y n < •)}(«) - £{Pr(K n+1 < ■)}(«) 



(e eA se« + A)" +1 ' 

The proof is thus complete. □ 

Lemma 6. Let m be an positive integer. For any x > 0, when e — > 0, E [/3J 1 ] — >■ 
E[7V r] . 

Proof. Since there is almost surely a finite number of points in [0, x], for almost all 
sample-paths, there exists r\ > such that AX, > rj for any j = 1, • ■ ■ , iV^. Hence 
for e < rj, Po(x) — N x . This implies that /?o(aO tends almost surely to N x as e goes 
to 0. Moreover, it is immediate by the very definition of (3q(x) that /3q(x) < N x . 
Since for any m, E [N™] is finite, the proof follows by dominated convergence. □ 



n > 0. 
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Let Li t (z), z,t s R, z < 1, be the polylogarithm function with parameter t, 
defined by 

OO 

fc=l 

For m a positive integer, consider the function of x 

OO 

(4) M™ g : x i > E [^(z)] = £ ^-(i). 

i=0 

Its Laplace transform is given by: 

p OO 

C{Mf Q }(s)= / E[/3 (o;r]e- si ^. 
Jo 

Corollary 7. Lei a 6e defined as follows: 

a = — — se . 
A 

TTie Laplace transform of the m-th moment of /3q(L) is: 

(5) C{M^}{s) = -^—LL 1 ' 



s(a + l) \a + l, 
which converges, provided that a > 0. 

Proof. Applying the Laplace transform of both sides of Eq. ((3]), we get: 

OO 

£{M™}( S )=^z m £{ K }( S ) 



OL T ■ ( 1 

-Li_„ 



' '^se es + l 



s(a+l) m V« + 1 , 
concluding the proof. □ 

We define |™| as the Stirling number of second kind [13) . i.e., |t| i s the 
number of ways to partition a set of m objects into k groups. They are intimately 
related to polylogarithm by the following identity (see |14| ) valid for any positive 
integer m, 

(6) u-«w=E — (1 _ z) l fc+1 



fe=0 

Corollary 8. TTie m-t/i moment of the number of clusters on the interval [0, L] is 
given by: 

(7) M%(L) = E {T} (7 - ( Aee ~ £A ) fe 
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Proof. Using §§§ in the result of Corollary [7J we get: 

m (-ir+^!{'T+ 1 }(i+^ 

l 0oS v ) ^ a k+i ( a + l) 

-Iv^ -1 

/ , Cfe m i, i 

s ' — ' or 

where the coefficients Ck, m are integers given by: 



771 /• 

Using the following identity of Stirling numbers |15| , 



wc,m — 

j=k 



m + 1 I / j 
J + l 



j=0 ^ 



o. 



we find that co, m = for m a positive integer. So we can write the Laplace transform 
of the moments as 

/c=i 

and apply the inverse of the Laplace transform in both size of Eq. ([3]) to obtain: 



£^i(L-fc £ )*(Ae--)% 



fc=i 



>M 



According to Lemma when e — ► 0, we obtain 



M%(L) = E [iV™] = £ ^L(LA) fc l {L>0} . 
fc=i 



Hence, for any A > 0, 

re / t- X Nfc-, I W | /r -vlfc-i 

>-{L>0}, 



fc=l ' fe=l ^ J 



which shows that 



fc! 



Thus, we have proved ((JJ for any positive integer m. □ 
Theorem 9. For any n, L, A and e, we have: 

{L/ej-n 



(8) Pr(A,(L) = n) = i E - (n + i)e)Ae^) 

i=0 
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Proof. Since (3q{L) < Nl and since E [e sAfL ] is finite for any s £ R, we have, for 
any s > 0: 



E 



e sf3 {L) 



OO 

£(-l) fe ^E[/3*(L)] 

k=0 



Rearranging the terms of the right-side hand and substituting MJ^(L), by the result 
of Eq. , we obtain: 



E 



(s) j { 3 



fe=0 \ j=k J 




Furthermore, it is known (see [T5]) that 



j=k J >• > 



Hence, 



E 



fc! 



e-*W| = ]T(L-fce) fe (A e 



L {L>ke}- 

k=Q 



By inverting the Laplace transforms, we get: 

" (-i)V^ 



where S a is the Dirac measure at point a. After some simple algebra, we find the 
expression of the probability that an interval contains n complete clusters: 

Pr(ML) =n)= p n (L) = ^ ~( n + W^T*, 

i=0 

concluding the proof. □ 

Lemma 10. For x > 0, p n (x) has the three following properties: 

i) p n {x) is dijferentiable; 
ii) lim x ^ 00 p„(a;) = 0; 
Hi) limx^oc dp 2i x) = 0. 

Proof. Let j be a non-negative integer. The function is obviously differentiable 
when x/e ^ j. Besides, we have 



lim p n {x) — lim p n (x) = lim — - — ( (x — (n + j)e) — 

X — ►ej - * - x— >-ej~ x — >e^' + J. 



n+j 



Since the right-hand term function of x is zero as well as its derivative for all j, 
the function is also derivable when x/e = j, which proves i). Items ii) and hi) are 
direct consequences of Final Value theorem in the Laplace transform of p n and its 
derivative. □ 

The expression of p n gives us a Laplace pair between the x and s domains: 

o) te> v t±r {x - {n+i)e) l) <u — ae +1 . 

V ' n! i\ V <V {ase" + l) n+1 

We can use this relation to find the distributions of Bi and U n . 
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Theorem 11. The distributions of Bi andU n , respectively fsi{x) and fu n {%) are 

d 



Xe-^poix - e) + e eA —p (x~e) 



l {x>e}, 



(10) 
and 

(11) fu n ( x ) = Ae~ eA p„_i(a; - e)l {x>e} , 

where the expressions of po(x — e) and -£;Po(x — e) are straightforwardly obtained 
from Eq. ©. 

Proof. According to Corollary [3J 

w[o - sBi i _ 1 (A + s) 



A ^-se SiL + 1 



= Ae 



-e\' 



A s± se ^ + l 



= Ae- eA e- £S £ {p (-)} (s) + e^e^sC {#,(•)} (s). 

Here, using the inverse Laplace transform established in Eq. and remembering 
that po(x~) = 0, we get an analytical expression for fBi(x), proving Eq. (fTU|) . 

Proceeding in a similar fashion, we can find the distribution of U n by inverting 
its Laplace transform given by Corollary 0] 

1 



E[e 



-sU n 



( s^L 



se se + 1 



Ae" 



se es + 1 



We thus have Eq. ([TT]l. 



□ 



We can also obtain the probability that the segment [0, L] is completely covered 
by the sensors. To do this, we remember that the first point (if there is one) is 
capable to cover the interval [X% — e, X\ + e]. 

Theorem 12. Let R m ,n{x) be defined as follows: 



b/ej-l 



( e -\ey+ n ^ (A[(l -i)e-x\y 
3=0 J ' 



Then, 



(12) Pr([0, L] is covered) = Rq,i(L) - e- Xe R 0A {L - e) 

- e- Ae i?!, (L) + e- 2A£ i?i,o(i - e). 
Proof. The condition of total coverage is the same as 

{Vze [0, L], 3Xi 6 [0,L] [X x - e,Xi + e]}, 
which means that: 

{[0, L] is covered} {B x > L - X t } n {ATi < e}. 

Hence, 

Pr([0, i] is covered) = / Pr(Bi > L — X\\X\ = x)dP Xl (x), 
Jo 
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and since B\ and X\ are independent: 

re />oo 

Pr([0, L] is covered) = / / f Bl (u)\e~ xX dudx. 
Jo J L-x 

The result then follows from Lemma [10] and some tedious but straightforward al- 
gebra. □ 

4. Other Scenarios 

The method can be used to calculate p n for other definitions of the number of 
clusters. We consider two other definitions: the number of incomplete clusters and 
the number of clusters in a circle. 

4.1. Number of incomplete clusters. The major difference with Sec.[3]is that a 
cluster is now taken into account as soon as one of the point of the cluster is inside 
the interval [0, L\. So, for instance, in Fig. O we count actually n + 1 incomplete 
clusters. We define f3' (L) as the number of incomplete clusters on an interval [0, L\. 

Theorem 13. Let G(k) be defined as 

G(k) = (-if (V- £ m_BL _ 1{T>ke} 

for keN+ and G(-l) = e~ AL . Then 

Pr(#,(L) = n) = J2 (-iy +n r)(G(i-l) + G(t)), forn>0. 



Proof. The condition of /3q(L) > n is now given by: 

f {AX + U n -x <L) ifn>l, 
[{AX <oo} if ra = 0. 



{P'o > n} 
We define Y n as 

Y n = 



AXo + Un-! ifn>l 
if 7i = 0. 



Repeating the same calculations, we find the Laplace transform of Pr(/3g(.) = n): 

A e eA e es 



£{Pr(^(.)=n)}( S )= S + A " (^-e» + l) 



- if n > 1, 



' if n = 0. 



A + s 

With this expression, following the lines of Lemma O we obtain 

— r rn + 11 1SI 1 A / Ae AeX 



r -i- 1 "i 

j — 1 L J 



s & A + s 



fe=i 

Then, we write: 

A 1 _ (-l) fe 1 k 1 7 1N k ~ l 



ill s i \ \ 



X + ss k A^ 1 A + s ^s 1 \ A 

2 — 1 V 

to find an expression with a well known Laplace transform inverse, and after in- 
verting it, we obtain: 



fc=0 ' 
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Expanding the Laplace transform of the distribution of /3' (L) in a Taylor series and 
rearranging terms, we get 



E 



sj3' (L) 



1 + G(0)£ 



r j 



/ i + 1 



,fe=i 



Now, we use another recurrence that Stirling numbers obey |15) . 

-(fe + 1) 



i + i 
fe + i 



j 

k + i 



to get: 



E 



x- 7 j j + 1 



IT 



1 



+ (*+!) 



(e »-i)* + _(e--l) 



J 

k + l 

k+l 



k\ 



k\ 



Hence, 



E 



sf)' (L) 



= l + J2(G(k-l) + G(k))(e- s -l) k 



k=l 



Inverting this expression for any non-negative integer n, we have the searched dis- 
tribution. □ 

4.2. Number of clusters in a circle. We investigate now the case where the 
points of the process are deployed over a circumference and we want to count the 
number of complete clusters, which corresponds to calculate the Euler's Character- 
istic of the total coverage, so we call this quantity \- Without loss of generality , 
we can choose an arbitrary point to be the origin. 

Theorem 14. The distribution of the Euler's Characteristic, x{L), when the points 
are deployed over a circumference of length L is given by 



(13) Pr( X (L) = n) = e-^l {n=0} + (1 - e~ XL ) 



eX YL/e\-n 



E 

i=0 



(-1)* 



([L-(n + i)e]Ae- £A ) 



e\\n-\-i— 1 



L + (n + i) l - - < 



1 



{X(L) > n} & 



for n > 0. 

Proof. If there is no points on the circle, x{L) = 0. Otherwise, if there is at least 
one point, we choose the origin at this point and we have equivalence between the 
events: 

{[/„-i + B n <L}n{N L >0} ifn>l, 
{AX <oo} if n = 0. 

In Fig. ||]we present an example of this equivalence. 
We can define Y n as 

U n -i + B n if n > 1 
if n = 0, 

to find the Laplace transform or Pr(x(L) = n): 
(14) r{Pr( X (0 = n)}(s) = (1 - e^ L ) A + s 



Yn = 



A A 



(^e» + l) 
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Bi 



B 2 



B r 



B 



n+1 



t 



Dr 



N L 




B r , 



Figure 4. Illustration of the condition equivalent to x(X) > n. 
Since the coverage of the last point on [0, L] overlaps the cluster 
with a point in zero, they are actually contained in the same cluster 



The number of clusters is almost surely equal to the number of points when e — > 
so 



E [ X (Lr\ = (1 - e~ XL )\e-^ £ \™\([L- ke]X 

fe=i L ^ ' 



eA\fe-l 



\e- £A ) 



L + k[ --€) ) l{L>fc £ } 
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Expanding the Laplace transform in a Taylor series and rearranging terms, as we 
did previously, yields 



E 



-«X(i) 



(l-e- XL )\e- eX Y^ 



k=0 



(\L - ke]\e 



k-l 



L + k 



1 



L {L>ke} 



E 

j=k 



(sy r j 



Since 



2^ j\ \k 



j=k 



k\ 



we can directly invert this Laplace transform, add the case where there are no 
points for x(£) = 0, and the theorem is proved. □ 



5. Examples 

We consider some examples to illustrate the results of the paper. Here, the 
behavior of the mean and the variance of fto(L) as well as Pr(/3o(L) — n) are 
presented. 

From Eq. Q, we have that E [/3o(L)] is given by: 

E [#,(£)] = (L-6)Ae- £A l {L>e} . 

This expression agrees with the intuition that there are three typical regions given a 
fixed e. When A is much smaller than 1/e, the number of clusters is approximatively 
the number of sensors, since the connections with few sensors will unlikely happen, 
which can be seen from the fact that E [/3q(L)] — > LX when A — > 0. As we increase 
A, the mean number of direct connections overcomes the mean number of sensors 
and, at some value of A, we expect that E [/3q(L)] decreases, when adding a point 
is likely to connect disconnected clusters. We remark that the maximum occurs 
exactly for e = 1/A, i.e., when the mean distance between two sensors equals the 
threshold distance for them to be connected. At this maximum, E [f3o(L)] takes the 
value of (L/e — l)e _1 . Finally, when A is too large, all sensors tend to be connected 
and there is only one cluster which even goes beyond L, so there are no complete 
clusters into the interval [0, L]. This is trivial when we make A — > oo in the last 
equation. Figure \5\ shows this behavior when L = 4 and e = 1. 
The variance can be obtained also by Eq. ([7]): 

Var(A)(£)) = (L- e)Ae- eA l {i>£} + (L — 2e)X 2 e - 2iX l {L>2e} 

IT \2\2 -2eA-i 

- [L - e) A e 1{l>c}, 

and under the condition that L > 2e: 

Var(/3 (L)) = (L - e)\ e - eX + e(3e - 2 J L)A 2 e~ 2e \ 

Fig. [5] shows a plot of Var(/3o(£)) in function of A for L — 4 and e = 1. We can 
expect that, when A is small compared to e, the plot should be approximatively 
linear, since there would not be too much connections in the network and the 
variance of the number of clusters should be close to the variance of the number of 
sensors given by XL. Since /3o(L) tends almost surely to when A goes to infinity, 
Var(/3o(i)) should also tend to in this case. Those two properties are observed in 
the plot. Besides, we find the critical points of this function, and again, A = l/e is 
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E[/3q] in function of A 




012345678 



A 

FIGURE 5. Variation of the mean number of clusters in function 
of A when L = 4 and e = 1. 

c§ Var(/?o) in function of A 




01234 5 678 

A 

FIGURE 6. Behavior of the variance of the number of clusters in 
function of A when L = 4 and e = 1. 

one of them and at this value Var(/3o(£)) = (L/e)e~ 1 + (3 — 2L/e)e~ 1 . The other 
two are the ones satisfying the transcendant equation: 

xp L — e 

Ae = 

2e(2L - 3e) 



By using the second derivative, we realize that l/e is actually a minimum. Besides, 
if L < 2e, there is just one critical point, a maximum, at A = l/e. 
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The last example in the section is performed with the result obtained in Theo- 
rem IH1 We consider again L = 4 and e = 1 to obtain the following distributions: 

Pr(A)(i) = 0) = l-3Ae- A + 2A 2 e- 2A -l/6A 3 e- 3 \ 

Pr(0b(L) = l) = 3Ae- A -4A 2 e- 2A + l/2A 3 e- 3 \ 

Pr(A)(L) = 2) = 2A 2 e~ 2A - l/2A 3 e- 3A , 

Pr(A)(i) = 3) = l/6A 3 e- 3A , 

Pr(/3 (i)>3) = 0. 

Those expressions are simple and they have at most four terms, since L = 4e. We 
plot these functions in Fig. [7] The critical points on those plots at A = 1/e are 
confirmed for the fact that, in function of A, for every n, Pr(x(L) = n) can be 
represented as a sum 

i=0 

where the coemcients qij are constant in relation to A. However, (Ae~ Ae )* has a 
critical point at A = 1/e for all i > 0, so this should be also a critical point of 
Pr(x(L) = n). If A is small, we should expect that Pr(x(L) = 0) is close to one, 
since it is likely to N have no points. For this reason, in this region, Pr(x(L) = n) 
for n > is small. When A is large, we expect to have very large clusters, likely 
to be larger than L, so it is unlikely to have a complete cluster in the interval and, 
again, Pr(%(£) = 0) approaches to the unity, while Pr(x(L) = n) for n > become 
again small. 




Figure 7. Probabilities of connectiveness, Pr(/3o(£) = n), for n = 
0, 1, 2, 3, in function of A when L = 4 and e = 1. 
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